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ABSTRACT 

ITie (Joddard A|xillo-Soyuz Goody n;iniic-s Kxjx;rimenl was performed to demon- 
strate the feasibility of traekinn and recovering high frequencj comix>nenls of 
the earni's gravity field by utilizing a synchronous orlnling iracking station such 
as A'I'S-(j. It was in reality the culmination of an effort started at Goddard in 
19G7/G8 to utilize synchronous orbiting tracking stations for NASA missions. 
(Vonbun and Mengei ]‘JG7). Gravity anomalies of say ■'» MGl.S or burger having 
wavelengths of 300 to 1000 kilometers on the earth's surface are im|x)rtant for 
geologic studies of the upper layers of the earth's crust. Using a low orbiting 
U- 230 km) spacecraft, such as A|x)llo-.S)yuz continuously tracked liy a synchro- 
nous satellite such as ATS-G, has for the first time detected from space short 
wavelength Earth's gravity anonuilies (Vonbun 1071). Two prime areas of data 
collection have been selected for this particular exiierimcnt. One area is the 
center of the Africjui continent and the second area is the Indian Ocean Deiires- 
sion centered at 5' north latitude and 7.G'' east longitude. 

Preliminary results show that the detectability objective of the exi)erimvnt lias 
been met in Ixith areas as well as at several additional anomcUous areas around 
the globe. Gravity anomalies of the Karakoram and Himalayan mountain ranges, 
ocean ti’enches, as well as the Diamantina Depth, cxiuld "he seen" to quote some 
specific examples. K\irther analyses are planned to demonstrate the utilization 
of the range rate signatures for the actual recovery of gravitational anomalies 
and their distribution on the earth's surface. 
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IN’ ri{()i)Uc:Ti()N 

'ITie objCL’tivfS of the A|Xillo/Soyuz (ieodymuiiies KxiJeriment were: 

(a) To demonstrate the detectability of short waveleiif^h (i.e. .‘IdO km and 
larjrer) features of the Karth's gravity field, 

(b) To evaluate the "IliKh/Low" satellite to satellite tracking (SST) concept 
for {feodynamics applications 

anil 

(c) To test the recoverability of short wavelength features of the earth's 
gravity field. 


The A|X)llo-Soyuz flying in a low orbit as mentioned, is particularly sensitive to 
gravity luiomalies. For example, the ratio of sensitivity to the 2.'Sth degree terms 
of the earth's gravity ixitential having a wavelength of say ItiOO kilometers is 
approximately 13 to 1 when going from a (IKOS-S t>n>e orl)it (810 km) to the low-er 
(230 km) A|x)llo orbit. ITius this mission gave us ;ui excellent opj)ortunity, which 
we could fortunately fully ex])loit, to test the satellite-to-satellite tracking 
concept. 

NoiTnaJ satellite perturbation methods used in the past for (he computation of 
the earth's gravity field cannot lie used for such short v/avelengths because this 
would require a spherical harmonic exyiansion to an order :uid degree 120 or 
larger. 


1. MAl’lIKMATICAL MKI’llOUS 

In general, the earth's gravitational field is in most cases represented in terms 
of a finite spherical harmonic series expansion whose coefficients have been de- 
termined from a combination of sati Hite ground tracking data and surtacc gravi- 
metric data. As is known, such a limited representation of the g.obal geoix)tential 
does not adequately "model" local variations in the gravity field, 'fhe reason 
being is that these local variations are too small in wavelengths and therefore 
would requir.' a field exj;ansion to a ve- y high degree an order, i.e., in the order 
of 100 to 200 which seems impractical 'f not imfxjssible. On the other hand, 
ralTer small local gravity anomalies are imix)rtant for geologic investigations of 
the upper crust as well as for studies of the ccean toixjgraphy (Vonbun lO?.")). The 
ex])eriment as described and executed gives 's an opjxjrtunity to improve our 
knowledge of the Karth's gravitational field jn 'ocalizcd areas by "directly" 
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measuring locaJ gravity variations by detecting small velocity changes of the low 
orbiting spacecraft. "Small" in this sense means 1 to 10 millimeters j)er second. 
Before going into some details on the gravity anomaly analyses some basics about 
the tracking system will lx; presented. Figure 1 shows the schematic of the ATS-6/ 
Apollo links. As can Ije seen, the Apollo is connected wiih Madrid via two links, 
namely, the ATS-6 link and a direct ground link. 'ITiis provides a total loop as 
far as the measurements are concerned which will be 'jsed in the future to leduce 
the errors in the system but will not be discussed in this particular report. Here 
we shall concentrate only on the link Mad -id-ATS-G-Aixjllo and back. Figure 2 
depicts those quantities which are actually measured during the experiment. 'Hie 
system is constructed in such a fashion to measure the total range rate lietween 
the ground station, in this case Madrid, and the Apollo spacecraft in earth orbit. 

The total range rate is the sum of the range rates between Madrid-ATS (rj ), 

ATS- Apollo (r^ ), Aix)llo-ATS (r^) and ATS-Madrid (i^) (Schmid, Vonbun 1974). 

For our analysis we will, however, only concentrate on the range rate variations 
between Apollo and A'PS that is ^ which is just a dot product of the spacecraft 
velocity vector and the unit vector between Apollo and A're as seen in f'igure 2. 


( 1 ) 

One of the first tasks during the experiment was to evaluate the noise character- 
istics of the total Lystem. Analysis done prior to the ex|)eriment predicted a 
range rate noise level of approximately 0. 05 cm sec utilizing a 10 second integra- 
tion time depicted in Figure 3. This figure further shows that the actual experi- 
mental range noise level was lower than the predicted one is in the order of 
0.03cm ^sec. These values were obtained using a small ixjrtion of an orbit so 
that possible orbit modelling errors where minimized. It should further Ije 
noted that these noise values are rather consistent through all the ex|)erime.ital 
arcs that have been analyzed. 

In this paper, two approaches have been taken to test the detectability of gravity 
anomalies. The first one was to assume a single mass anomaly and to compute 
in a simple analytical fashion, the expected variation of the spacecraft velocity 
due to this anomaly, (Vonbun 1972). 'fhe velocity variations in the radial and 
horizontal components are then: 
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see Fig. 4. 

'I'he above equations assume that the spacecraft travels directly over the anomaly 
and are only valid up to ±5° in true anomaly which corresponds to about 1.5 
minutes of Qight lime, 'fhe dashed curve on Figure 4 shows the radial velocity 
variation due to a 5 mgal anomaly (4° by 4°). The size of the anomaly block was 
roughly twice the orbital height and represents the area of sensitivity of a space- 
craft in orbit to an anoma'y on the ground, (Schwartz 1970). Because of the 
extreme low noise in the data and the real geometry (spacecraft is not directly 
over the extended mass rather than point mass), a computer simulation was 
necessary to improve on the alx)ve shown analytical results. 

'fhe dotted line in Figure 4 depicts the result of this computerized analysis, 
which will now b? briefly described. As .shown in Figure 2 and mentioned earlier, 
the value is equd (Bryan 1974): 


f Fj + Tj 4 4 
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this is the range rate as actually measured at the Madrid ATSR station. We are, 
however, only interested in the range rate between the A{X)llo spacecraft and the 
ATS spacecraft. That is: 
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which is the clot product of the spacecraft orhitjil velocity vector v’ and the unit 
vector tjetween A|»llo and ATS-G. Since the ATS-6 orbit is almost independent 
of any gravity variations because of its extren.c height (3G, ()()() km) its range rate 
can simply be subtracted from the measured one. ITiat is 

^ ■ - / (ATS Madrid) 

which brings us to one quantity needed for our further analyses, niis range rate 
is on the other hand a function of 

/a KX,. Xa,,. Cl) (G) 

where Xa and Xa-p^ are Ajxillo and AT’S-G state vectors respectively and Cl 
represents the gravity field, TTie major aim in this exp- •'iment was to find the 
relationship between surface gravity anomalies g and the Aixj11o/AT’S-G range 
rate variations To be more specific, we are interested only in the varia- 

tion in ,'a* thus we want to vary equation (G). 
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and also: 


(7) 


= [( Avpy + Avp ) ^ = (ohsL-rval-cominilod ) K.miiv K.itv 

The sensitivity cjoefficient ( IXa/ ‘(j) is now contained in the equation 7 together 
with A Cl. Helating the actual measurements with the gravity jx)tential variation 
A G. Since we are observing local anomalies over very short time only as can be 
seen in Figure 2 say in the order of a few minutes the first two terms on the 
right side of the series exjiansion can be neglected thus equation 7 reduces to 
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IMfferentiatinK (9) with resjx'ct to CJ yields: 



G G G 


( 10 ) 


Since ( X^/i*G) is our variable to be determined we rewrite equation (10) in the 
following manner assuming intorchangealiility of differential operators; that is: 


{ (I^ h / rV \ / ■'V \ 
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( 11 ) 


WTiere V' = U + T and U represents the normal p<)tential of the earth's gravity 
field in spheiical harmonics and 


T 



AgS(,) d.’ 


Where Ag = AG + Ag^^^ using Ag^^^ as the local average anomaly ol a block 
(say 4° X 4° as an example), is the disturbing ixjtential due to the anomalies 
A g, U is the radius of a sphere whose volume is equal to the volume of the spheroid 
and S(vj) the Stokes IXmction. 

Introducing U and T into (11) yields finally: 
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Kquation (12) is now solved for the sensitivity coefficient (d X^/d G) by numeri- 
cal integration. Using equation (8), the above computed sensitivity coefficient is 
correlated and the relationship between \i^ and AG established. In actuality 
the following procedures were used in detecting the gravity anomalies. 

An orbit of 1 to 1.5 revolutions was computed for the Apollo spacecraft utilizing 
the (Preliminary Goddard Solution) PGS-108 gravity field using all tracking 
available, that Is ground station as well as SST tracking data. From these Apollo 
orbits together with the ATS-6 orbit, the range rate between ATS-6 nnd Apollo 
were computed. These range rate data were subtracted from the actual measured 
ones resulting in theAp^ as shown in Figures 6 and 8. The reason that these 
values of A/j^ are finite rather than zero is the fact that we do not know the 


gravity fiehi to a very high accuracy which wcjuld necessarily include all the 
gravity lUiomalies. I'his is of course the essence of this particular ex|>eriment. 


II. PHKLIMn^ARY ItKSULTS 

Kxj)eri mental data were collected in the Eastern liemispheric region where the 
A|)ollo-Soyuz was in range of the A'I’S-6 satellite which is iwsitioned over the 
ecjuator at approximately .'15” east longitude to perform the Indian Communica- 
tions Exjjeriment. The duration of each A}X)llo orbit, visible to ATS, was approx- 
imately 50 minutes or slightly greater than one half ot an orbital revolution. 'Hie 
actual SST data passes used were only about 40 minutes long in order to assure 
non-disturbance of these data due to the earth's atmosphere. All data were suc- 
cessfully collected using the A'I'S-C spacecraft and the Madrid Al'S- h station as 
the prime ground station. Preliminary results confirm that the range rate noise 
of 0.05 cm/s computed Ijefore the exfieriment was exceeded by the experimental 
data which show an average noise of approximately 0.03 cm/sec (Fig. 3). Fig- 
ure 5 shows the ground track of revolution 115 of the .\tx)llo which coincidentally 
passes over two anomalies in the African area, as indicated. The gravity varia- 
tions shown resulted from differencing two Goddard Earth Models (GEM-7 and 
GEM-G). The even GEM are (;onstructed by using ground and satellite data where 
the odd numbered models are constnicted only from satellite data. 'Phis differ- 
ence thus reflects in essence the surface gravity anomalies though somewhat 
filtered by the models. 'ITiese two anomalies have a variation of about ±5 mgals. 
In Figure 6 are depicted the theoretical variations of a 5 mgal anomaly (center! 
and sui)erimiK)sed (lower part) resulting in a total characteristic signature one 
would expect. In addition, the figure on the right also depicts the observed 
radial velocity variations. As can be seen, a fairly good agreement exists be- 
tween the simulated and observetl variations of the radial component of the space- 
craft velocity except for a time delay. The dashed line shows finally the time 
shifter! simulated variation of Ap., 

For this analysis the PGS lOH Field of order and degree of 25 was used. In ihat 
particular case the gravity anomalies are in the onier of 5 mgals and extending 
over 300 to 2000 km were actually observed since the wavelength of the gim\1ty 
anomaly and the resolution of this field are roughly compatible. This shows the 
usefulness of satellite-to-satellite tracking data to resolve anomalies of say 
t5 mgals on the Earth's surface. 


The second test area, namely the Indian Ocean Depression, is shown in Figure 7 
in some detail utilizing again the Goddard P^arth Morlel, Gem 0, but augmented 
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with 1° by 1° surface gravity data n Figure 8, four orbital passes, namely 8, 23, 
83, ana 113, are plotted In the upper part of the graph. These orbits pass as 
can be seen over the Indian Ocean Depression as well as over the Himalayan 
Mountain range. The lower portion depicts the actual measurement residuals 
obtained, that is, the observed minus computed values. Both, the Indian Ocean 
Depression and anomalies in the Himalayan area can easily be "seen" in the actual 
data. Taking the Indian Depression as an example, the shape of the range rate 
varlatlon.s do actually fit the expected one computed in a similar fashion as those 
shown in Figure 0. Detectaljillty of the gravity field perturbations have thus been 
demonstrated, manifesting themselves in the expected form of small spacecraft 
velocity variations. The peak values of these variations could in most cases be 
closely correlated with anomalies shown on the detailed map of the global geoid. 
Other features may have been detected but cannot be verified at this tim»» since 
there is a lack of surface gravimetric data in these areas. This is particularly 
true for the southern portion of the globe. One may however assume that the 
detected "gravity signatures" in the southern part of the earth rep^'esent existing 
anomalies not presently known too well. 

Orbit determination required for the proper analysis of the experimental data to 
date involves the use of the satellite-to- satellite range rate data together with 
the direct ATS-6 tracking ?md the Unified S-Band tracking data. Error analysis 
erformed (Indicates that one and two revolution arcs seem to be the best in 
order to minimize the effeids of unmodelled error in the detection and the re- 
covery of gravity anomalies. Longer arcs seem to introduce larger errors because 
the influence of the Earth's gravity field errors onto the orbit errors become more 
and more pronounced as the arc length increases. Shorter arcs however have 
also rather large errors since not enough tracking data are available. ITius, 
the use of one to two orbit revolutions for the analyses seems to be a good opti- 
mum. Air drag corrections have been applied for the orbit determination of the 
Apollo- Soyuz. 'Fhe satellite-to-satellite tracking data were only used during 
those periods which were judged free of on-board propulsion activities. Since 
the antenna is about 3 meters separated from the spacecraft axis any attitude 
motion (.-xjrrection) of the spacecraft, that is a rotation of the spacecraft, trans- 
forms itself into a range rate between the Apollo-Soyuz and the ATS-6. One 
could therefore quite wrongly interpret a spacecraft rotation as a gravity anomaly 
on the ground. Detailed analysis of attitude control has not been performed at 
this time. Attitude data are, however, available and will be u.sed in the future to 
further refine the range rate residuals as measured betv'een Apollo-Soyuz and 
ATS-6. 
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CONCLUSION 

Basetl upon the forcKolng, It can be stated that the " .nh-low" satelltte-to- 
satellite tracking technique can be utilized for earth gravity anomaly detection. 
Anomalies whose dimensions are eqtial or longer than the height of the low 
spacecraft and have values in the order of 5 mgals have been detected with a 
good 10:1 signal to noise ratio. The range rate of 0.05 cm/s conservatively com- 
puted prior to the experiment was actually exceeded by the measured noise of 
0.03 cm/s. It was possible for all the passes checked so far, that the measured 
range rate residuals could be correlated with known geoidal features of short 
wave length of say 300 km or larger. It was further demonstrated that the range 
rate signatures could be repeated orbit by orbit thus Increasing the confidence 
in the experimental data. 
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Figure 8. Itange liate Signatures of ASTP Revolutions in the 
Indian Ocean and Himalayan Areas. 
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